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ABSTRACT 


The  time  equation  of  satellite  orbit  theory  is  defined  as  the  equation 
used  to  predict  the  time  at  which  a  satellite  will  either  cross  the  earth's 
equator  or  pass  through  the  perigee  point  of  its  orbit,  depending,  respec¬ 
tively,  on  whether  the  nodal  or  anomalistic  period  is  used.  This  paper 
gives  a  derivation  of  the  time  equation  through  terms  of  order  of  the 
fourth  time  derivative  of  the  period.  It  is  shown  that  when  the  time  equa¬ 
tion  is  expressed  as  a  polynomial  the  initial  period  is  the  sum  of  the  co¬ 
efficients  of  the  variable.  Since  the  time  equation  is  also  used  to  deter¬ 
mine  the  numerical  values  of  its  coefficients  by  least  square  fitting  the 
equation  to  observations,  equations  are  derived  for  the  first  and  second 
time  derivatives  in  terms  of  these  coefficients.  An  analysis  of  some 
NORAD  and  NSSCC  data  indicates,  at  least  for  certain  satellites,  that 
determining  optimal  values  of  these  coefficients  to  use  in  the  satellite 
position  predicting  method  is  closely  related  to  the  time  interval  over 
which  the  least  square  method  is  applied. 
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SECTION  1  -  INTRODUCTION 


The  time  equation  of  satellite  orbit  theory  is  defined  as  the  equation 
used  for  predicting  the  time  at  which  a  satellite  will  either  cross  the 
earth's  equator  or  pass  through  the  perigee  point  of  its  orbit,  depending, 
respectively,  on  whether  the  nodal  or  anomalistic  period  is  used.  This 
paper  gives  a  derivation  of  the  time  equation  through  terms  of  order  of 
the  fourth  time  derivative  of  the  period.  By  the  use  of  Kepler's  third 
law  the  time  derivatives  of  the  period  can  be  replaced  by  time  derivatives 
of  the  semi-major  axis,  but  this  is  not  carried  out  in  the  paper.  It  is 
shown  that  when  the  time  equation  is  expressed  as  a  polynomial  the  initial 
period  is  the  sum  of  the  coefficients  of  the  variable.  Since  these  coef¬ 
ficients  are  usually  obtained  by  least  square  fitting  the  time  equation  to 
observations,  the  period  and  also  the  first  and  second  time  derivative  of 
the  period  are  expressed  in  terms  of  these  coefficients. 


SECTION  2  -  THE  TIME  EQUATION 

The  time  equation  can  be  written  in  the  forms  (1)^,  (2),  or  (3). 

Tn  =  Tq  +  Pq  ^1  ^2  ■*■  •  •  •  %-l  • 

Tn  =  To  +  2  A.  •  F.  (N).  (2) 

j=l  •'  ■' 

Tn  =  To  +  Bj  •  N-j.  (3) 

We  will  call  (1)  the  exact  time  equation,  where  is  the  time 

required  for  a  satellite  to  complete  its  Nth  revolution  since  the  epoch 
time  Tq.  If  the  periods  Pj  (j  =  0,  1,2,  N-1)  are  the  nodal  periods, 

then  Tj^  is  the  time  of  the  mh  (south  to  north)  equatorial  crossing.  If  Fj 
is  the  anomalistic  period  then  T^  is  the  time  at  which  the  satellite  passes 
through  perigee. 

figures  in  parentheses  refer  to  equations. 
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Eqiiation  (2),  while  not  in  general  use,  is  given  in  this  p>aper  because 
it  serves  as  an  internnediate  equation  in  deriving  (3)  from  (1),  and  also  be- 

por 

cause  the  coefficients  Aj  represent  a  single  L  ^  defined  in  equation  (4) 
following.  In  (2),  Aj  is  the  epoch  period  Pq. 


Equation  (i)  is  most  generally  used  to  obtain  It  has  been  used 

for  several  years  at  the  National  Space  Surveillance  Control  Center 
(NSSCC),  now  the  Space  Track  Research  and  Development  Facility,  L.  G. 
Hanscom  Field,  Bedford,  Massachusetts,  and  it  has  also  been  used  at  the 
North  American  Air  Defense  Command  (NORAD),  Colorado  Springs, 
Colorado.  Both  NORAD  and  NSSCC  have  found  that  fairly  accurate  ephem- 
erides,  for  most  artificial  earth  satellites,  can  be  issued  when  either  the 
quadratic  or  cubic  equation  in  N  is  used.  The  coefficients  B;  in  (3)  are 

pqr 

expressible  as  a  sum  of  the  L  qiiantities  defined  in  (4). 


In  Appendix  I  the  Taylor  series  for  is  used  in  conjunction  with 
(1)  to  obtain  the  following  expressions  for  Aj  .  Fj  (N).  To  simplify  the 

pqr 

writing  of  these  expressions  we  first  define  L  by  the  equation 

0 1> 


^_pqr  ^  P  .  , 
St  0  0 


*  •  ♦ 

Po 

r 

«  •  •  * 

Po 

s 

•  •  •  • 

Po 

21 

31 

41 

(4) 


where  the  dots  denote  time  derivatives.  We  see  from  (4)  that  t  is  either 
zero  or  one,  and  for  the  latter  value  it  follows  q  =  r  =  s  =  0.  The  terms 
of  (2)  for  j  =  1,2,  3,  ...,  12  are  given  in  (5)  through  (16). 


Fi(N)-Ai  =  N  L^o°0°- 

(5) 

N  (N-l) 

F2(N)-A2  *  -  • 

2 

(6) 

N  (N-l)  (2N-1)  L^oO 

F3(N)-A3  = 

6 

(7) 

_  ^  N  (N-l)  (N-2)  L  0^0 

^^(Nj-A^  =  ^ 

(8) 

2  2  300 

N  (N-l)  L  iQ 

F.(N)-A  =  -  . 

5  ^  A 

(9) 

2 


APGC-TDR-62-8 


F^(N).A6  * 
F^{N)•A^  * 
F8(N)-A8  = 
F9(N)-Ag  = 
Fio<n)-Aio 


Fii(N)-Aii  = 


N  (N-1)  {N-2)  (N-3)  lVo° 


24 


N(N-l)  (N-2)  (2N-1)  L^oo 


N(N-l)  (2N-1)  (3N^  -  3N-1)  l'^qi 


30 

N  (N-1)  (N-2)  (N-3)  (N-4)  L 


140 

00 


120 


N  (N-1)  (N-2)  (42N^  -  39N-1)  l\^q 
120 

N  (N-1)  (N-2)  (22N^  -  69N  +  29)  L^qq 
120 


Fi2(N).Ai2  * 


N  (N-1)  (N-2)  (16N^  -  22N  +  2)  L^qq^ 
120 


(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


In  the  L^^*'notation  the  highest  power  of  N  multiplying  a  given  1-'^^^*' 

isp  +  q  +  r  +  s  +  t.  When  the  sums  of  these  superscripts  and  subscripts 
are  the  same,  as  in  (12)  through  (16),  then  numerical  values  given  by 
these  equations  are  under  normal  circumstances  of  about  the  same  order 
of  magnitude.  Hence,  if  we  wish  to  retain  only  terms  through  the  third 
time  derivative,  then  (12)  through  (16)  can  be  omitted.  If  we  retain  only 
terms  through  order  of  the  second  time  derivative,  then  (9)  through  (11) 
can  also  be  omitted. 


The  coefficients  in  (3)  are  obtained  by  collecting  like  powers  of  N  in 
(5)  through  (16),  and  are  given  in  (17)  through  (21). 

110  201  120  130  211  400 

_  .  100  ^00  ^00  ^00  L,  00  L-  00  ^01 

^  on  ■■  - - - - - -  - 

^  2  6  3  4  3  30 


140 
'  00 


310 
^  10 

"60~ 


29L 


221 

00 


60 


302 
-  00 

30 


(17) 
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B 


2 


B 


3 


B4 


no  201  120  300  ,,  130  211 

L-OO  L.00  L,  00  L.  10  llL  00  00 

“1  2  2“  ^  4  ^  24  6 


140  310  221  302 

5L  00  5L  10  15L  00  5L  qo 


12  8  8  12 


201  120  300  130  _  211  400 

^00  ^00  ^10  ^00  '^00  ^01 

“3“  ^  “6  2  4  6~  ^  “3 

140  310  221  302 

7L  00  5L  10  7L  00  5L  00 

+  -  +  -  +  -  +  - , 

24  3  3  6 


300  130  211  400  140  310 

L-  10  ^00  ^00  ^01  ^00  *^^10 

4  24  3  2  12  8 


9  221  7  302 

■  ?  ■  ll  “O  • 

400  140  310  221  302 

^01  ^00  10  11^00  2L’qq 

Be  =  -  +  -  +  - +  -  +  - . 

^  5  120  20  60  15 


(18) 


(19) 


(20) 


(21) 


Equations  (17)  through  (21)  give  the  coefficients  in  (3).  Assume  that 
numerical  values  of  these  coefficients  have  been  obtained  for  a  satellite 
by  least  square  fitting  (3)  to  the  observations.  We  first  determine  the 
period  at  epoch  Tg  by  inspection  of  (17)  through  (21)  and  find  that,  if 
terms  of  order  of  the  s-1  time  derivative  of  the  periods  are  retained, 
we  have 


s 

Pg  =  .2  B;  (s  =  2,  3,  4,  5). 

j=l  J 


(22) 


We  next  determine  the  rate  of  change  of  the  period  Pg  in  terms  of  these 
coefficients.  Inspection  of  (17)  through  (21)  reveals  that  expressions  for 
Pq  become  complicated  quite  rapidly  as  s  increases.  For  s  equal  to  4 
or  5,  it  is  necessary  to  solve  a  cubic  or  biquadratic  equation,  respec¬ 
tively,  in  order  to  obtain  Pg.  Since  in  most  applications  of  (3),  as  at 
NORAD  and  NSSCC,  s  is  set  equal  to  2  or  3,  we  shall  derive  only  expres¬ 
sions  for  these  values  of  s.  For  s  =  2,  we  obtain  from  (17)  and  (18): 
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PqPo 


PqPo 


B2  . 


and  we  obtain  from  (23)  and  (24): 

Pq  =  Bi  +  B2. 


PO  = 


2B2  2B2 


+  ^2  Pq 


For  8  =  3,  we  have  from  (17),  (18),  and  (19): 

PqPo  PqPo^  Pq^Pq 


Bl  =  Pq 


12 


B2  = 


PqPo  PqPo^  Pq^p'o 


Bi  = 


PqPo^  Po^Pq 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 


From  (27)  through  (29)  we  obtain  Pq,  Pq,  Pq,  and,  since  these  quantities 
should  be  obtained  in  that  order,  the  following  equations  are  given  both  in 
terms  of  the  B's  and  of  the  previously  calculated  quantities. 

For  s  =  3,  we  have: 


Pq  =  Bi  +  B2  +  B3, 

'Bj  -  3B2  -  5B^ 


(30) 


Pq  =  1 


=  1  - 


Bi  +  B2  +  B3J 


1/2 


4B2  +  6B3 


nl/2 


1  - 


PO 


(31) 
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IOB3  +  2B2  -  2Bj  2(Bj  -  3B2  - 
fBj  +  B2  +  B3)2  *  (Bi  +  B2  +  B3)3/2 
6B3  Po^ 


(32) 


In  general,  sufficient  accuracy  of  Pq  is  obtained  if  (31)  is  expanded 
into  an  infinite  series  and  only  the  first  term  is  retained. 

2B2  +  3B3 

P  ,,  -  .  (33) 

PO 


SECTION  3  -  THE  TIME  EQUATION  AS  USED  AT  NORAD  AND  NSSCC 


The  method  of  predicting  satellite  positions  which  was,  until  recently, 
used  at  NORAD  and  NSSCC  was  originally  formulated  by  NSSCC.  In  this 
paper  the  NORAD/NSSCC  method,  or  simply  method,  refers  to  the 
method  used  until  recently.  Since  the  beginning  of  1962,  improvements 
have  been  introduced  at  SPADATS  which  include  a  differential  correction 
sequence  in  the  predicting  method  and  which  makes  the  system  semi¬ 
automatic  in  its  use. 

In  April  1961,  the  author  obtained  from  NSSCC  the  IBM  orbital  ele¬ 
ment  cards  for  all  satellites  in  orbit  from  January  1959  to  April  1961. 

The  elements  given  on  these  cards  are  nodal  elements.  These  IBM  cards 
contain  the  orbital  elements  used  to  obtain  the  Space  Track  bulletins 
(ephemerides)  issued  during  this  time  interval.  The  bulletins  were  gen¬ 
erally  issued  to  cover  a  time  interval  of  3  to  10  days.  One  of  the  pur¬ 
poses  for  obtaining  these  cards  was  to  undertake  an  error  analysis  in 
the  hope  that  this  analysis  would  shed  some  light  on  possible  ways  of 
improving  predictions.  While  improvements  appear  to  be  definitely 
possible,  at  least  for  certain  satellites,  the  author  believes  that  a  com¬ 
plete  proof  of  such  improvements  must  be  verified  with  the  observational 
data  at  NORAD  or  NSSCC.  Therefore,  the  author  has  made  arrange¬ 
ments  to  secure  the  necessary  observational  data  and  reduction  program 
from  NSSCC. 
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If  improvements  are  to  be  made  in  the  NORAD/NSSCC  method,  one 
should  first  determine  what  appears  to  be  the  major  cause  of  correctable 
errors,  and  then  determine  if  a  partial  or  complete  elimination  of  this 
cause  leads  to  better  predictions.  Before  pointing  out  what  appears  to  be 
the  major  cause  of  errors  in  predictions,  some  opinions  on  other  causes 
of  errors  will  be  given.  First,  the  author  believes  the  observational  data 
received  at  NORAD  and  NSSCC  are  sufficiently  accurate  to  permit  more 
precise  predictions.  Secondly,  the  equations  employed  in  the  NORAD/ 
NSSCC  method  are  sufficiently  accurate  to  permit  more  precise  predic¬ 
tions,  and  it  is  therefore  not  necessary  to  include  more  theoretical  terms 
in  these  equations  until  such  times  as  errors  in  predictions  are  explainable 
by  their  omission.  Thirdly,  while  unpredictable  variations  in  air  drag 
are  certainly  a  serious  drawback  to  long  range  predictions,  or  some¬ 
times  even  to  short  range  predictions  when* rapid  acceleration  variations 
occur  due  to  geomagnetic  disturbances  (such  as  the  October  and  November 
I960  events  [1]^),  it  is  possible  to  make  longer  range  predictions  and 
still  remain  within  the  limit  of  permissible  errors. 


The  author  believes,  and  knows  that  a  good  many  people  in  NSSCC 
and  NORAD  concur,  that  the  major  causes  of  errors  in  predicting  satel¬ 
lite  positions  are  errors  in  the  coefficients  of  the  time  equations.  In 
order  to  point  out  the  importance  of  these  coefficients,  the  first  three 
equations  used  in  the  NORAD/NSSCC  method  are  listed  below.  Since  in 
practice  it  is  generally  necessary  to  retain  the  cubic  term  in  the  time 
equation  only  during  the  decay  phase  of  a  satellite,  we  shall  therefore 
drop  this  term.  We  employ  the  notation  used  at  NORAD  and  NSSCC. 


Tn  =  To  +  (N  -  Nq)  Pq  +  (N  -  No)2 


UN  *  <^0  ^'^N  "  “0  +  (Tn  -  Tq) 


Cq. 

(34) 

2  ^ 

• 

(35) 

2 

2  UQ 

(36) 

• 

2 


The  only  assumptions  in  these  three  equations  are  that  the  period, 
the  right  ascension  of  the  ascending  node  n,  and  the  argument  of  perigee 
u  can  be  adequately  represented  by  their  Taylor  series.  This  assumption 
on  the  period  is  used  in  the  appendix  to  derive  the  time  equation.  In  (34), 
Pq  is  called  the  nodal  period  at  revolution  number  Nq,  and  Cq  = 


PqPo 


While  it  is  true,  as  pointed  out  in  the  previous  section,  that 


2 


Numbers  in  brackets  refer  to  references. 
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when  the  quadratic  equation  is  used  the  actual  nodal  period  at  revolution 
Nq  should  be  assigned  the  value  of  Pq  +  Cq  ,  the  assignment  of  this  value 
will  contribute  nothing  to  the  accuracy  of  predictions  unless  the  error  in 
Pq  is  less  than  the  magnitude  of  Cq.  In  general  this  is  not  the  case. 
However,  if  (34)  is  fitted  to  observations  and  the  Pq  and  Cq  are  obtained 
to  such  accuracy  that  the  magnitude  of  Cq  affects  the  accuracy  of  the 
initial  period,  then  the  nodal  period  should  be  assigned  the  value  Pq  4  Cq, 
and  this  value  used  in  computing  the  anamolistic  period,  the  semi-major 
axis,  etc.  Of  course,  values  of  Pq  obtained  by  the  least  square  fitting 
of  (34)  to  the  observations  are  to  be  retained  in  (34),  and  we  must  never 
replace  Pq  by  Pq  +  Cq  in  the  time  equation.  In  what  follows  we  shall 
always  refer  to  Pq  in  (34)  as  the  initial  period  at  revolution  number  Nq. 

In  predicting  future  positions  of  a  satellite,  (34)  is  first  used  to  ob¬ 
tain  Tjij,  and  then  T|^  -  Tq  is  substituted  into  (35)  and  (36)  to  obtain  Ajyj 
and  Hence,  it  is  readily  seen  that  any  errors  in  the  time  equation 

produce  errors  in  subsequent  computations  in  the  predicting  method.  It 
is  therefore  essential  that  the  time  equation  be  made  to  predict  as  accu^ 
rately  as  possible. 

In  order  to  understand  how  errors  are  introduced  into  the  coefficients 
of  the  time  equation,  we  shall  point  out  how  the  coefficients  of  the  time 
equation  are  usually  obtained.  Let  us  first  rewrite  the  time  equation  in 
the  form, 


T  (N  +  Ni)  =  T  (Ni)  +  N  P  (Ni)  +  C  (Ni).  (37) 

Assume  that  at  epoch  revolution  the  T,  P,  and  C  on  the  right  hand 
side  of  (37)  are  known.  Also  assume  that  we  restrict  the  time  residuals 
AT  (N  -f  N^)  (differences  between  observed  and  predicted  equatorial 
crossings)  to  a  maximum  absolute  value  of  30  seconds.  The  dots  in  Fig.  1 
represent  the  time  residuals  AT.  After  N  revolutions  (N  =  Nj  -  N^)  it 
becomes  apparent  that,  unless  corrections  are  made  in  the  time  equation, 
the  value  of  AT  will  exceed  30  seconds  by  revolution  number  Nj^.  Hence, 
it  is  imperative  that  a  new  Space  Track  bulletin  be  issued.  The  least 
square  method  is  then  used  to  fit  the  time  equation  to  the  observations, 
let  us  say  over  the  time  interval  T(Nj)  -  T(N^),  and  a  new  set  of  coeffi¬ 
cients  obtained  for  (37).  If  the  epoch  revolution  Nj  is  selected  for  the 
next  bulletin,  then  to  conform  to  the  notation  used  in  (37),  must  be 
replaced  by  Nj  in  (37).  Inherent  in  the  application  of  the  least  square 
method,  as  applied  to  obtain  the  new  time  equation  at  Nj,  is  the  expecta¬ 
tion  that  the  time  residuals  will  continue  to  lie  along  the  solid  curve  in 
Fig.  1  if  the  predictions  are  allowed  to  continue.  If  the  time  residuals 
did  continue  to  lie  along  this  curve,  then  we  should  expect  accurate  pre¬ 
dictions  when  the  coefficients  at  Nj  are  used  in  the  time  equation. 
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Fig.  1:  Time  Residuals  Versus  Revolution  Number 


However,  in  many  cases  the  time  residuals  do  not  continue  to  lie  along 
this  curve,  and  in  these  cases  the  least  square  method  does  not  give  us 
a  set  of  optimal  coefficients.  The  time  residuals  could,  in  fact,  continue 
along  the  dashed  curve  and  actually  reach  a  negative  maximum  permis¬ 
sible  error  before  reaching  a  positive  maximum  permissible  error.  A 
serious  dilemma  is  then  encountered  in  the  problem  of  making  more  ac¬ 
curate  predictions,  since  we  are  faced  with  the  almost  impossible  task 
of  determining  along  which  curve  the  time  residuals  will  lie. 

The  author  is  currently  investigating  a  number  of  methods  which, 
for  certain  satellites,  give  more  precise  predictions.  In  order  to  de¬ 
scribe  one  of  these  methods,  suppose  that  the  time  equation  (37)  has  been 
least  square  fitted  to  the  incoming  observations  over  the  time  interval 
T(Nj)  -  T(Nj),  and  we  then  use  (37)  with  Nj  replaced  by  Nj  to  predict 
future  nodal  crossings.  Why  was  the  time  interval  T(Nj)  -  T(N^)  used  to 
obtain  the  orbital  elements  at  T(Nj)  rather  than  some  other  time  interval? 
The  answer,  of  course,  is  that  the  orbital  elements  at  T(Ni)  were  used 
over  this  time  interval  to  obtain  the  time  residuals,  and  we  have  no  other 
time  interval  to  use.  This  answer  leads  to  the  following  question.  Is  it 
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possible,  for  certain  satellites,  that  the  time  equation  should  be  least 
square  fitted  over  a  much  longer  time  interval  in  order  to  obtain  more 
accurate  coefficients  and  hence  minimize  the  time  residuals?  The  author 
does  not  have  available,  at  the  present  time,  the  necessary  data  to  com¬ 
pletely  answer  this  question.  However,  it  is  possible  to  approximate  the 
coefficients  which  would  be  obtained  by  least  square  fitting  over  a  longer 
time  interval  by  using  the  epoch  revolutions  and  epoch  times  on  three 
successive  sets  of  orbital  element  cards.  As  an  example  of  this  proce¬ 
dure  we  use  Space  Track  Object  No.  63,  I960  Pi  1,  Tiros  II.  From  the 
three  successive  sets  of  orbital  element  cards  issued  at  epoch  revolu¬ 
tions  560,  650,  and  910,  we  obtain  the  respective  epoch  times  T(560), 
T(650),  and  T(910).  These  data,  in  conjunction  with  the  time  equation 


T  (N  +  Nj)  =  T  (Nj)  +  N  P  (Nj)  +  C  (Nj),  (38) 

give  the  two  equations 

T(560)  =  T(910)  -  350  P(910)  +  (350)^  C(910)  and  (39) 

T(650)  =  T(910)  -  260  P(910)  +  (260)^  C(910),  (40) 

which  can  be  solved  for  P(910)  and  C(910).  The  equation 

T(N  +  910)  =  T(910)  +  N  P(910)  +  C(910)  (41) 

is  then  used  to  predict  the  epoch  times  on  later  sets  of  orbital  element 
cards.  Before  using  (41)  to  predict  these  later  epochs,  we  set  N  =  -  470 


and  predict  the  epoch  T(440).  The  difference  between  this  predicted 
epoch  and  the  epoch  on  the  orbital  element  cards  is  three  seconds,  and 
this  indicates  that  the  values  of  P  and  C  in  (41)  should  permit  fairly 
precise  predictions  of  future  epochs.  The  results  of  such  predictions 
are  given  in  Fig.  2,  where  the  AT  curve  gives  the  time  residuals  which 
would  be  obtained  if  the  P  and  C  given  on  the  orbital  element  cards  issued 
at  revolution  number  910  were  used  in  (41),  and  the  AT  curve  gives  the 
time  residuals  which  would  be  obtained  if  (39)  and  (40)  are  solved  for  P 
and  C  and  these  values  used  in  (41). 

The  decision  as  to  whether  a  longer  time  interval  should  be  used  for 
a  specific  satellite  should  certainly  not  rest  on  a  single  example  as  shown 
in  Fig.  2.  This  decision  should  be  made  on  the  consistency  of  improve¬ 
ments.  For  Tiros  II,  when  the  method  was  applied  to  14  consecutive 
sets  of  cards,  including  the  sets  issued  at  epoch  revolution  910,  in  each 
application  we  had  AT  <AT.  For  another  example  we  use  Space  Track 
Object  No.  29,  I960  Beta  2,  Tiros  I.  In  Fig.  3  the  sets  of  cards  issued 
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at  epoch  revolution  numbers  100,  250,  and  490  were  used  to  obtain  the 
coefficients  of  the  time  equation  at  490,  In  Fig.  4,  the  sets  of  cards 
issued  at  revolution  numbers  2500,  2650,  and  2970  were  used  to  obtain 
the  coefficients  of  the  time  equation  at  2970.  In  Figs.  3  and  4,  the  AT 
and  AT  notation  is  the  same  as  in  Fig.  2. 

A  comparison  of  the  AT  and  AT  curves  in  Figs.  2  through  4  indicates 
that,  for  certain  satellites,  the  accuracy  and  also  the  time  intervals  over 
which  Space  Track  bulletins  can  be  extended  are  dependent  on  how  much 
of  a  satellite's  recent  history  is  incorporated  into  determining  the  co¬ 
efficients  of  the  time  equation. 

In  order  to  further  point  out  that  it  is  the  coefficients  of  the  time 
equation  which  are  critical  in  accurate  position  predictions,  we  use  data 
from  Space  Track  Object  No,  28,  I960  Beta  1,  rocket  body  of  Tiros  I. 
Starting  at  epoch  revolution  250,  the  equation 

•  • 

7  “2  50 

‘^N  =  “250  +  -  ”^250)  “250  +  -  T250)  - .  (42) 

2 

with  T250  *  109.  69590;  01250  =  185.  7047;  W250=  4.  117376; 

“250  . A 

— ~ — s  9»  34  X  10"  ,  predicts  that  values  of  oij^  on  the  successive  29  sets 

of  orbital  element  cards  with  a  maximum  error  less  than  0,  025  degree. 

In  making  these  predictions  the  values  of  Tj^^  on  the  respective  cards 
were  used.  These  cards  cover  a  time  interval  of  351  days,  or  5100 
revolutions,  and  during  this  time  the  perigee  has  moved  through  an  angle 
of  1446.  5  degrees.  The  accuracy  of  (36),  when  applied  to  I960  Beta  1, 
is  not  to  be  construed  as  indicative  of  the  accuracy  with  which  (36)  will 
predict  the  argument  of  perigee  for  other  satellites  over  long  time  inter¬ 
vals.  For  example,  when  (36)  is  least  square  fitted  to  1958  Alpha  data 
(see  [2])  over  a  time  interval  of  approximately  two  years,  the  maximum 
(aj  residuals  are  approximately  two  degrees.  The  example  cited  here, 

1958  Beta  1,  is  solely  for  the  purpose  of  pointing  out  that  while  other 
equations  employed  in  the  NORAD/NSSCC  method  are  sufficiently  accu¬ 
rate  for  longer  range  predictions,  this  fact  is  of  little  value  unless  one 
can  predict  with  sufficient  accuracy  the  value  of  Tjy|.  In  obtaining  the 
AT  curve  in  Fig,  5  the  coefficients  of  the  time  equation  at  epoch  revolu¬ 
tion  490  were  obtained  by  using  the  epoch  times  given  on  successive  sets 
of  orbital  element  cards  issued  at  revolution  numbers  100,  250,  and  490. 
The  AT  curve  was  obtained  by  using  the  P  and  C  given  on  the  cards  is¬ 
sued  at  epoch  revolution  490. 

Inspection  of  the  AT  curve  reveals  that  the  maximum  error  in  pre¬ 
dicting  4860  revolutions  into  the  future  is  probably  about  240  seconds. 
While  this  maximum  error  is  approximately  eight  times  our  assumed 
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permissible  time  error,  the  use  of  these  predicted  values  in  (42), 
instead  of  the  values  on  the  orbital  element  cards,  will  only  change 
the  predicted  values  by  approximately  0.01  degree.  This  shows  that 
whereas  errors  in  time  may  be  excessive,  such  errors  do  not  necessarily 
introduce  serious  errors  into  the  angles. 

It  should  be  strongly  pointed  out  that  the  method  of  obtaining  the  co¬ 
efficients  of  the  time  equation,  in  order  to  obtain  the  AT  values  given  in 
Figs.  2  through  5,  is  not  recommended  by  the  author  for  operational  use, 
because  the  coefficients  obtained  in  this  manner  are  generally  not  the 
optimal  values  (see  Fig.  6).  What  is  recommended,  however,  is  that, 
in  order  to  make  more  accurate  and  longer  range  predictions,  serious 
consideration  be  given  to  the  time  interval  over  which  the  least  square 
method  is  applied. 

The  AT  curve  in  Fig.  6  for  Space  Track  Object  No.  36,  I960  Epsilon 
3  (metal  object  from  Sputnik  IV),  was  obtained  by  using  in  the  time  equa¬ 
tion  the  P  and  C  on  the  orbital  element  cards  issued  at  epoch  revolution 
number  4440.  The  AT  curve  was  obtained  when  the  three  successive  sets 
of  orbital  element  cards  issued  at  revolution  numbers  4150,  4280,  and 
4440  were  used  to  calculate  the  P  and  C  at  revolution  4440.  While  AT 
<AT  in  Fig.  6,  the  ST  curve  is  certainly  not  the  optimal  time  residual 
curve.  If  the  period  P  used  in  obtaining  the  AT  curve  is  changes  by  0.  25 
second,  the  AT  curve  is  obtained.  The  AT  curve  is  approximately  a  sine 
curve  with  a  period  of  27  days  and  an  amplitude  of  12  seconds. 

The  St  curve  can  be  considered  as  an  empirical  curve  which  is 
drawn  after  the  time  residuals  are  known.  However,  if  the  criteria 
(referred  to  in  the  following  section)  are  applied  to  the  P  and  C  values 
used  in  the  AT  curve,  the  cubic  term  in  the  time  equation  will  only  account 
for  approximately  five  seconds  of  the  error  at  the  end  of  the  AT  curve. 
Hence,  since  we  are  justified  in  using  the  quadratic  time  equation  over 
this  time  interval,  the  assumption  that  the  residuals  are  caused  by  errors 
in  the  P  and  C  values  and  by  the  omission  of  periodic  terms  leads  to  cor¬ 
recting  the  P  by  0.  25  second,  and  gives  the  best  fit  curve  which  is  the  AT 
curve. 


SECTION  4  -  COMMENTS 


In  Section  2  and  Appendix  I,  the  time  equation  was  derived  to  include 
all  secular  terms  of  order  of  the  fourth  time  derivative  of  the  period.  The 
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Fig.  6;  Satellite  I960  Epsilon  3  -  Metal  Object  from  Sputnik  IV 


previous  sections  could  leave  the  reader  with  the  impression  that  it  is 
only  P  and  P  that  are  important  in  the  time  equation,  and  this  is  certainly 
not  the  case.  The  terms  containing  P  through’p  'in  (17)  through  (21)  are 
important  in  making  an  accurate  error  analysis  of  applications  of  the  time 
equation.  Due  to  the  fact  that  accurate  expressions  for  P  through  P  can 
be  derived  in  terms  of  P  and  P,  the  author  has  developed  criteria  for  de¬ 
termining  good  approximations  to  the  optimal  time  interval  over  which  the 
time  equation  should  be  least  square  fitted.  Applications  of  these  equations 
have  increased  the  accuracy  of  predicting  satellite  positions  and  have  also 
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improved  satellite  decay  predictions.  In  addition  to  P  and  P,  it  may 
sometimes  be  necessary  to  include  the  eccentricity  and  its  first  time 
derivative  in  these  equations.  The  author  is  continuing  to  develop  and 
evaluate  these  equations  and  expects  to  publish  the  results  at  a  later  date. 

In  deriving  the  time  equation,  periodic  terms  were  intentionally 
omitted  because  they  were  not  included  in  the  NORAD/NSSCC  method. 

Their  omission  is  entirely  justifiable  because  the  residuals  are  primarily 
due  to  inaccuracies  of  the  secular  terms  and  to  variations  in  atmospheric: 
density.  Actually,  the  residuals,  once  the  bad  errors  in  secular  terms 
are  eliminated,  do  show  up  in  nearly  all  the  curves  given  in  Section  3. 

The  decision,  therefore,  on  whether  to  include  or  not  include  the  periodic 
terms  depends  primarily  on  the  accuracy  to  which  the  secular  terms  have 
been  obtained.  For  most  orbiting  satellites  the  inclusion  of  periodic  terms 
in  the  NORAD/NSSCC  method  should  not  be  an  exceedingly  difficult  task. 
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APPENDIX  I 

DERIVATION  OF  THE  TIME  EQUATION 


In  this  appendix  the  time  equation  is  derived  to  include  all  terms  of 
order  of  the  fourth  time  derivative  of  the  period.  The  derivation  applies 
both  to  nodal  anc  anomalistic  periods  as  defined  previously.  The  period 
P^  is  expressible  as  a  continuous  function  of  time  Tj,^. 

Pn  =  P(Tn).  (1-1) 

and  the  Taylor  series  for  the  period  is 

Pn  =  PQ  +  -  Tq]  Pq  +  [Tn  - 

•  •  •  • 

+  [Tj(J  -  Tq]^-^  +  [Tj^  -  Tq]  - +  ....  (1-2) 

31  41 


The  exact  time  equation  is 

Tn  =  To  +  Pq  +  Pi  +  Pz  +  •••  +  Pn-1,  (^-3) 

where  Pn_i  is  the  time  required  for  a  satellite  to  complete  the  Nth  revol¬ 
ution.  From  (1-3)  we  have 


Tl  -  Tq  =  Pq* 


(1-4) 


Substitution  of  (1-4)  into  (1-2)  gives  Pj, 
Pi  into  (1-3)  gives 

Z” 
PO  PO 


Tz  -  Tq  =  2Po  +  PqPo  + 
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and  substitution  of  this  value  of 
‘V  •  •  4*  •  •  • 

.  ,,.5, 

3'  4' 


Equation  (1-5)  is  now  substituted  into  (1-2)  to  obtain  It  is  necessary 

to  repeat  the  above  procedure  until  Fq  is  obtained.  Since  the  calculations 
become  lengthy,  the  results  are  given  in  the  following  tabulated  form, 
where 


(1-6) 
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The  numbers  in  the  jth  row  are  the  coefficients  of  the  quantities 

for  the  Pi  in  that  row.  The  expressions  for  the  coefficients  of  the  LP*!*" 

8 1 

quantities  for  Pp^-l  easily  obtainable  as  functions  of  N  but  are  unnec¬ 
essary  for  our  purpose.  Inspection  of  the  table  shows  that  the  sum  of  the 


100 

00 


numbers  in  the  second  column  is  N,  and  therefore  the  coefficient  of  L 
in  the  time  equation  is  N,  In  the  third  column  we  have  1  +  2  +  3  +  . , .  + 

In  the  fourth 


(N-1),  and  therefore  the  coefficient  of  i®  N_(N_^) 

00  2 

column  we  have  1^  +  2^  +  3^  +  +  (N-1)^,  hence  the  coefficient  of  L-^qq^ 

N(N-l)  (2N-1)  ,  ,  3  3  3  3 

is  - ,  In  the  fifth  column  we  have  1+2^  +  3'^  +  ,,,  +  (N-1)  , 


and  therefore  the  coefficient  of  i® 


300..,  (N-1)^ 


.  In  the  sixth  column 


we  have  1^  +  2^+  3^  +  .,.  +  (N-1)^,  and  the  coefficient  of  becomes 

N(N-l)  (2N-1)  (3N^  -  3N-1) 

30 

The  coefficients  of  the  remaining  seven  quantities  are  obtain¬ 

able  in  the  following  manner.  Since  it  was  only  necessary  to  obtain  Pq 
in  order  to  derive  the  coefficient  of  .  we  shall  illustrate  the  method 

by  obtaining  the  coefficient  of  Applying  the  method  of  differences 

to  column  9  of  Table  1,  we  have: 


1 


15 


35 


70 


126 


10 


20 


35 


56 


10 


15 


21 


1  1 

and  therefore  the  sum  of  N  numbers  in  the  ninth  column  is 

4N  (N-l)  6  N  (N-1)  (N-2)  4  N(N-l)  (N-2)  (N-3) 

^  21  3!  41 

N  (N-1)  (N-2)  (N-3)  (N-4) 


(1-7) 
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1 40 

Since  there  are  N-4  numbers  in  the  column  for  L'  qq  •  replace^N^in 
(1-7)  by  N-4,  When  this  substitution  is  made  the  coefficient  of  L  is 
found  to  be 

N  (N-1)  (N-2)  (N-3)  (N-4) 

Fq  (N)  s  - .  (1-8) 

9  120 

The  remaining  coefficients  are  found  in  a  similar  manner,  and  all 
12  of  the  coefficients  are  given  in  (5)  through  (16)  in  Section  2, 
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tivr»  tn  tBrm*  of  those  coofficionts.  An  onxlyiis  of  aomo  NORAD  and  tiv*#  in  terms  of  those  coefficients.  An  analysis  of  some  NORAD  and 

NSSCC  data  indicates,  at  least  for  certain  satellites,  that  determining  N85CC  data  indicateo,  at  least  for  certain  satellites,  that  determining 

optimal  values  of  these  coefficients  to  use  in  the  satellite  position  pre>  optimal  values  of  these  coefficients  to  use  in  the  satellite  position  pre¬ 
dicting  method  is  closely  related  to  the  time  interval  over  which  the  least  dieting  method  is  closely  related  to  the  time  interval  over  which  the  least 

square  method  is  applied.  square  method  is  applied. 


